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Abstract. We study the long Josephson junction (LJJ) model which takes into account the 

second harmonic of the Fourier expansion of Josephson current. The dependence of the static 
[~t | magnetic flux distributions on parameters of the model are investigated numerically. Stability 

CZ2 , of the static solutions is checked by the sign of the smallest eigenvalue of the associated Sturm- 

Liouville problem. New solutions which do not exist in the traditional model, have been found. 

Investigation of the influence of second harmonic on the stability of magnetic flux distributions 

for main solutions is performed. 
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1. Motivation and model 

The physical properties of the magnetic flux in Josephson junctions (JJ) are the contemporary 

superconducting electronics base. For the traditional JJ the current-phase dependence is a 

sinusoidal function. Such a model is described by the sine-Gordon equation. For a sufficiently 

► wide class of JJ the superconducting Josephson current as a function of magnetic flux ip can be 



t^J- ■ represented as a sine series [H [2] : 

oo 
Is = I c sin ip + \^ Im sin rrup . (1) 
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It was recognized recently that the higher harmonics in this expansion are important in many 
applications, in particular, in junctions like SNINS and SFIFS, where S is a superconductor, I 
is an insulator, N is a normal metal and F is a weak metallic ferromagnet [HE]- The interesting 
properties of LJJ with an arbitrarily strong amplitude of second harmonic in current phase 
relation were considered in J3j. 

Using only first two terms of the expansion (TJJ leads to the double sine-Gordon equation 
(2SG) [2j. 

— ip + a\ sin ip + a,2 sin 2tp — 7 = , x E (— I; I) . (2) 

Here and below the prime means a derivative with respect to the coordinate x. The magnitude 
7 is the external current, I is the semilength of the junction, a\ and 02 are the normalized 
amplitudes of the first and second harmonics of the Josephson current [U [5] . All the magnitudes 
are dimensionless. 



The boundary conditions for ([2]) have the form 

<p'(±l) = h e 



(3) 



where h e is external magnetic field. 

Numerical solution of the nonlinear boundary problem ([2]), (|3|) is solved on the basis of the 
continuous analog of Newton's method [7]. 

Stability and bifurcations of static solutions <p(x,p), where p = (l,ai,a2,h e ,y) are analyzed 
on the basis of numerical solution of the corresponding Sturm-Liouville problem [6] : 



ijj" + q(x)ip = Xip, ip (±Z) = 0, q(x) = aicostp + 2ci2Cos2(p. 



(4) 



The minimal eigenvalue \o(p) > corresponds to the stable solution. In case \o(p) < solution 
<p(x,p) is unstable. The case \o(p) = indicates the bifurcation with respect to one of the 
parameters p. We characterize the solutions of equation (J2J) by number of fluxons N(p) which 
is defined as 



N(p) 



1 
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tp(x) dx . 



(5) 



2. Deformation of the Meissner solution Mi 
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Figure 1. The internal magnetic field of the Meissner solution Mq for LJJ with 2/ = 10, 7 = 0, 
h e = 2 and a\ = 1 at different values of the parameter 02 . Left figure shows <p ' at negative a<i , 
right one shows at positive 02- 



In the "traditional" case ai = two trivial solutions (p = and <p = it of ([2~]), ([3]) are known 
at 7 = and h e = 0, which are denoted by Mo (N[Mq] = 0) and M v {N[M^\ = 1), respectively. 
Accounting of the second harmonic ai sin 2(p leads to the appearing of two additional solutions 
ip = ±arccos(— ai/202) denoted as M± ac (N[M± ac \ are not integer numbers and depend on the 
value of second harmonic). Stability properties of trivial solutions in dependence on parameter 
ci2 are considered in [7j and [8]. 

All solutions with N[ip] =0 we denote here by Mq, even they are changed by the influence of 
the external magnetic field and parameter 02- Within the LJJ model, the basic Meissner solution 
Mq demonstrates the screening of the external magnetic field. Deformation of Mq at different 
values of parameter 02 under the external magnetic field h e = 2 is shown in Fig. [TJ At a<i < 
(left panel) the screening effect diminishes when the coefficient a2 decreases and it amplifies at 
02 > (right panel) with growing of 02- 




Figure 2. The internal magnetic field of the distribution <J> at 21 = 10, 7 = 0, h e = 2 and 
a\ = 1 at different values of the parameter ai- Left figure shows </?' at negative 02, right one 
shows at positive 02- 



3. Deformation of the fluxon solutions $ n 

The solutions with iVjyj] = 1 which are not ip = 7T at h e = and 7 = we denote as $ 1 . 
The deformation of the (f'(x) of fluxon states under influence of the parameter a<i E [— 1; 1] was 
considered in [8j. We observe a qualitative change in the local minima with decrease of a<i in 
the interval [-1;0] and in the local maxima with increase of 02 in the interval [0;1]. With change 
of the coefficient 02 the number of fluxons corresponding to the distribution & 1 is conserved [6] 
i.e. dN/da,2 = and iVf^ 1 ] = 1. As external magnetic field h e is growing, more complicated 
stable fluxon with n = N[ip] = 2, 3, . . . appear which are denoted by <3? n . 

The effect of the second harmonic contribution on <1> 1 was studied in [3] and [9]. It was shown 
in [1] that taking into account at h e = the second harmonics in current phase relation with 
02 < —0.5 leads to the appearance of the "small" fluxon state additionally to the traditional 
"large" one. Taking into account (0), we call the solution at ai = —0.7 in Fig. [TJ left panel as 
Mq (in [1] it is a "small" fluxon), because iV[Mo] = 0. Number of fluxons of "large" fluxon is 
equal to N [large] = 1 and we denote it by $ . In [{J|, stability properties of "large" fluxon $ 1 at 
nonzero h e have been studied. Two coexisting stable 3? 1 like fluxons were demonstrated in some 
region of magnetic field. Here, we investigate this effect in case of two-fluxon and three-fluxon 
distributions. The relation between the "small" , "large" fluxons and trivial solutions under the 
influence of the external magnetic field and the second harmonic are the point of our further 
research. 

In Fig. [2] the internal magnetic field of the two-fluxon distribution $ 2 for 21 = 10, 7 = 0, 
h e = 2 and a\ = 1 at different values of the parameter 02 is presented. We observe a strong 
deformation of the fluxon distributions with decrease of the parameter 02 in the interval [-1;0]. 
At 02 = —0.5 the curve of internal magnetic field ^>'{x) has a plateau at the points x ~ —4, 
x ~ and x ~ 4 (Fig. [21 left panel). Further increase of the absolute value of 02 leads to 
the transformation of these plateaus to local maxima of the internal magnetic field. Thus, 
accounting of the 02 contribution qualitatively changes a shape of the fluxon distribution $ 2 . 
A similar deformation in the local maxima regions (points x ~ —2 and x ~ 2) is observed for 
02 > (Fig. [21 right panel). We stress that the number of fluxons is conserved iV[<I> 2 ] = 2 for 
all values 0,2- 



4. Stability analysis of the static fluxon distribution <l> n 

The two-fluxon solution <J? 2 in LJJ with 21 = 10 and 7 = is stable at 0.5 < h e < 2.5. 
Qualitatively, the behavior of the curves \o(h e ) for this solution is similar to the case of one- 
fluxon solution <1> 1 [9j. The only difference is that in the $ 2 case we don't have a region of 
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Figure 3. The dependence Xo(h e ) for <3? n at 21 = 10, 7 = 0, h e = 2 and ai = 1 at different 
values of the parameter a-i- Left figure shows Xo(h e ) for $ 2 at negative 02, right one shows 
Xo(h e ) for $ 3 at positive 02- 

coexistence of two stable branches that was observed for negative 02 in the & 1 case. 

Results for ai G [— 1;0] are presented in Fig. [3j left panel. When a^ decreases in (— 0.5; 0], 
the curve Ao(/i e ) moves to the right. At a^ < —0.5 the curve Xo(h e ) corresponding to the stable 
solution <£ 2 has two separate branches. 

The Xo(h e ) curves for 3> 3 distribution for 02 € [0; 1] are demonstrated in Fig. [31 right panel. 
It is seen, when 02 is growing in the interval [0; 1], the first bifurcation point moves to the left. 
Contrary, at a<i > 0.7 the bifurcation point moves to the right. The second bifurcation point 
moves to the left as a^ is growing from to ai w 0.7. Contrary, at 02 > 0.7 the second bifurcation 
point moves to the right. For three-fluxon state the following effect is observed. At 02 > 0.5 
new curves appear and we observe an existence of two different 3-fluxon states simultaneously. 

For trivial solutions our calculations show the reducing of Meissner screening when the second 
harmonic is negative. New fluxon solutions which appear in case 02 7^ have a probability to 
be observed in the experiment. So, it would be interesting to test this. 

As summary we note that our numerical investigations show that accounting of the second 
harmonic contributions significantly changes the shape and stability properties of trivial and 
fluxon static distributions in LJJ. 
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